A single Dirac-fermion Quantum Heat Engine 
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We studied the efficiency of two different schemes for a quantum heat engine, by considering a 
single Dirac fermion trapped in a one-dimensional potential well as the "working substance". The 
first scheme is a cycle, composed of two adiabatic and two iso-energetic reversible trajectories in 
configuration space. The trajectories are driven by a quasi-static deformation of the potential well 
due to an external applied force. The second scheme is a variant of the former, where iso-energetic 
trajectories are replaced by isothermal ones, along which the system is in contact with macroscopic 
thermostats. This second scheme constitutes a quantum analogue of the classical Carnot cycle. Our 
expressions, as obtained from the Dirac single-particle spectrum, converge in the non-relativistic 
limit to some of the existing results in the literature for the Schrodinger spectrum. 
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I. INTRODUCTION 

A classical heat engine consists of a cyclic sequence of 
reversible transformations over a "working substance", 
typically a macroscopic mass of fluid enclosed in a cylin- 
der with a mobile piston at one end 0, Hj]. The two 
most famous examples are the Otto and Carnot cycles. 
In particular, the classical Carnot cycle comprises four 
stages, two isothermal and two adiabatic (iso-entropic) 
ones. Ideal quasi-static and reversible conditions are 
achieved by assuming that an external force, which differs 
only inhnitcsimally from the force exerted by the inter- 
nal pressure of the fluid, is applied to the piston in order 
to let it move extremely slowly (HQ- On the other hand, 
the isothermal trajectories are performed by bringing the 
fluid contained by the cylinder into thermal equilibrium 
with external reservoirs at temperatures Tc < Th, re- 
spectively. 

A quantum analogue of a heat engine involves a 
sequence of transformations (trajectories) in Hilbert's 
space, where the "working substance" is of quantum 
mechanical nature [3T-flc|. One of the simplest concep- 
tual realizations of this idea is a system composed by 
a single-particle trapped in a one-dimensional potential 
well @] ■ The different trajectories are driven by a 
quasi-static deformation of the potential well, due to the 
application of an external force. Two different schemes 
of this process have been discussed in the literature, in 
the context of a non-relativistic particle whose energy 
eigenstates are determined by the Schrodinger spectrum 
I n this paper, we shall revisit these approaches, 
and we will study the performance of the correspond- 
ing heat engine for a single Dirac fermion. Since Dirac's 
equation describes the spectra of relativistic particles, 
results obtained for this case should reduce to the cor- 
responding ones from Schrddinger's equation in the non- 
relativistic limit. As we shall discuss below, the transi- 
tion between the relativistic and non-relativistic regimes 
is determined by the ratio X/L, with A = 2irh/mc the 
Compton wavelength of the particle, and L the width of 



the potential well. The non-relativistic limit corresponds 
to the regime where X/L <C 1, while evidence of the un- 
derlying relativistic nature of the spectrum manifests in 
terms of finite corrections in powers of X/L. Another 
limit of theoretical interest is the "ultra-relativistic" case 
of massless Dirac fermions, where X/L — > oo. An impor- 
tant realization of this former case in solid state systems 
is provided by conduction electrons in the vicinity of the 
so called Dirac point in graphene (Il1 - fl4| . 



II. A DIRAC FERMION TRAPPED IN A 
ONE-DIMENSIONAL INFINITE POTENTIAL 
WELL 

We define a one-dimensional potential well by the func- 
tion 



V{x) 



lim V [Q(x - L) + e(-x)}, 

Vo-*-oo 



(1) 



with 0(a;) the Heavisidc step function. The eigenvalue 
problem for a Dirac fermion trapped in this potential well 
is stated as follows 

Hrj> = Erj>, (2) 

where the Dirac Hamiltonian in 1+1 dimensions is [l5| 

H = -ihaid x +mc 2 a 3 + V(x)l. (3) 

Here, o\ and 63 are Pauli matrices, while ip{x) = 
(4>{x) , \{x)) is a two-component spinor. 

The solution, considering the boundary condition at 
the infinite barrier interfaces 4>(0) = 4>(L) = 0, is given 
by 

sin(n7ra;/L) 



M*) = A{_ iJ ^-- nizx/L) y (4) 

with associated energy eigenvalues organized in pairs 
with opposite signs [15| . 



E%(L) = ±mc' 



(y/l + (n\/2L) 2 -l) 



(5) 
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FIG. 1: (Color online)Pictorial description of the Force versus 
size of the potential well, for the E-Carnot Cycle described in 
the text. 



Here, we have substracted the rest energy, and A = 
2irh/ (mc) is the Compton wavelength. The positive (neg- 
ative) sign corresponds to the particle (anti-particle) so- 
lution, respectively [l5j . In what follows, we shall assume 
that a single particle is trapped by the potential well, and 
hence we shall keep the positive eigenvalue. 

The spectrum predicted by Eq.([S]) can be compared 
with the corresponding Schrodingcr problem, whose 
eigenvalues are given by 



E S n {L) 



n 2 TT 2 h 2 
2mL 2 ' 



(6) 



Here, a single eigenvalue for the energy is obtained and, 
moreover, it scales as n 2 , in contrast with the Dirac par- 
ticle case where a richer scaling with n is observed. In 
particular, in the regime A/L <C 1, we have: 



(n\/2L) 2 = E*(L), 



(7) 



which corresponds to the non-relativistic limit. Beyond 
this regime, relativistic corrections depending on the fi- 
nite ratio X/L are observed. Another interesting limit 
of Eq.© corresponds to a masslcss Dirac particle with 
A — > oo, where the spectrum reduces to the expression 



E n( L )\ m=0 



nithc 
L ' 



(8) 



This situation may be of interest in graphene systems, 
where conduction electrons in the vicinity of the so 
called Dirac point can be described as effective mass- 
less chiral particles, satisfying Dirac's equation in two 
dimensions Ill|-[l4l . 



III. A SINGLE-PARTICLE QUANTUM HEAT 
ENGINE 



As the "working substance" for a quantum heat engine, 
let us consider a statistical ensemble of copies of a single- 



FIG. 2: (Color online) Pictorial description of the Energy ver- 
sus Entropy (von Neumann) for the E-Carnot cycle. 



particle system, where each copy may be in any of the 
possible different energy eigenstates. We therefore say 
that the single-particle system is in a statistically mixed 
quantum state [16[ • The corresponding density matrix op- 
erator is p = t2 n p n (L)\il) n (L)}(il) n (L)\, with \tp n (L)} an 
cigenstate of the single-particle Hamiltonian Eq.Q, cor- 
responding to the spinors defined by Eq.((4]). This density 
matrix operator is stationary, since in the absence of an 
external perturbation [TJ| ihd t p = [H,p\ = 0. Here, the 
coefficient < p n {L) < 1 represents the probability for 
the system, within the statistical ensemble, to be in the 
particular state \ip n (L)). Therefore, the {p n (L)} satisfy 
the normalization condition 



Trp = ^2pn{L) = 1. 



(9) 



In the context of Quantum Statistical Mechanics, en- 
tropy is defined according to von Neumann [I?} as 
S = — fc^Trpln p. Since in the energy eigenbasis the den- 
sity matrix operator is diagonal, the entropy reduces to 
the explicit expression 



S(L) = -k B Y,Pn(L)ln( Pn (L)) 



(10) 



In our notation, we emphasize explicitly the dependence 
of the energy eigenstates {\ip n (L))}, as well as the prob- 
ability coefficients {p n (L)}, on the width of the potential 
well L. The ensemble- average energy of the quantum 
single-particle system is 

E=(H) = Tv(pH) = J2Pn(L)E n (L). (11) 



For the statistical ensemble just defined, we conceive 
two different schemes for a quantum analogue of a ther- 
modynamic heat engine. The first one, that we shall 
refer to as the E-Carnot cycle, consists on four stages 
of reversible trajectories: two iso-entropic and two iso- 
energetic ones. During the iso-energetic trajectories, 
the ensemble-average energy Eq. tjlip is conserved, while 
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during the iso-cntropic ones, the von Neumann entropy 
defined by Ea. ([TU|) remains constant. We distinguish 

this first scheme from the quantum Carnot cycle to ^ /'.'' ' a = 16 

be discussed next, where the iso-energetic trajectories 
in Hilbert's space are replaced by isothermal processes. 
During these stages, the system is brought into thermal 
equilibrium with macroscopic reservoirs at temperatures 
Tc < Th, respectively. 
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IV. THE E-CARNOT CYCLE 



0.2 > 



The system trajectories in Hilbert's space are assumed 
to be driven by reversible quasi-static processes, in which 
the walls of the potential well are deformed "very slowly" 
by an applied external force, such that the distance L is 
modified accordingly. Along these trajectories, the total 
change in the ensemble average energy of the system is 
given by 

dE = J2Pn( L )dE n (L) + J2dp„(L)E n (L) 

n n 

The first term in Eq. (|12[) represents the total energy 
change due to an iso-entropic process, whereas the second 
term represents a trajectory where the energy spectrum 
remains rigid. 

Let us consider first the iso-entropic process, where 
{p n (L)} = cnt. We remark that this represents a strong 
sufficient condition for the entropy to remain constant 
along the trajectory, but is not a necessary one[l8j]. Un- 
der quasi-static conditions, the external force driving the 
change in the width of the potential well is equal to 
the internal "pressure" of the one-dimensional system, 
F = —(dE/dL)s- Therefore, the work performed by the 
system against the external force, when the width of the 
potential well expands from L = L\ to L = L2, is given 
by 



dL 



8E\ 



9L J { Pn (L)} = {p n (L 1 )}=cnt. 

OO 

^2,p n (Lx) [E n (L 2 ) - E n (Li)\ . (13) 



n=l 



For the case of a Dirac fermion, the work performed 
under iso-entropic conditions is given after Eq. (|13[) and 
Eq.© by 



^ Pn (L 1 )(\/l + (n\/2L 2 y 



Yl + (nA/2i!) : 



(14) 



Notice that our sign convention is such that, for an ex- 
pansion process L 2 > L\, the work performed by the sys- 
tem is negative^, indicating that the ensemble-averaged 



FIG. 3: (Color online)Efficiency of the E-Carnot cycle, calcu- 
lated after Eq. (|28|l . as a function of the expansion parameter 
a = Lc / Lb- 



energy is decreasing during expansion, as in a classical 
ideal gas. 

Let us now consider an iso-energetic process, that is, a 
trajectory in Hilbcrt space defined by the equation dE = 
0. The solution to this equation, for L £ [L\, L 2 ], is given 
by the path 



J2Pn{L)E n {L) = ^p„(Li)£„(Li), 



(15) 



n=l 



along with the normalization condition Eq.(|9]). Clearly, 
by definition an iso-energetic process satisfies 



dE = SWi^2 + <5<2i->2 = 0, 



(16) 



with SW^ 2 = (SE) {pn{L]}=cnL and SQ^ = 
(SE)^ E fk)}=cnt ■ I n the former equation, the integral 
along the trajectory L\ — > L 2 gives 



AE = Wi^ 2 + Qi^2 = 0. 



(17) 



The first term W\^. 2 corresponds to the mechanical 
work performed by the system against the external force, 
which drives the change in the width of the potential well 
at constant energy. The second term Q\^ 2 = —W\^, 2 
corresponds to the amount of energy exchanged by the 
system with the environment, in order to rearrange its 
internal level occupation. This equation is in analogy 
with the first law of thermodynamics for macroscopic sys- 
tems, when considering a reversible process over a clas- 
sical ideal gas which is being compressed/expanded at 
constant internal energy conditions. The first term has 
a precise correspondence with the mechanical work for 
expansion/compression, whereas the second is in corre- 
spondence with the heat exchanged by the gas with the 
environment in order to satisfy total energy conservation. 

According to the previous analysis, the heat exchanged 
by the system with the environment along the iso- 



4 



energetic process is given by 



dL 



(18) 



Evidently, Eq. (jT5|) combined with the normalization con- 
dition Eq.© are not enough to uniquely define the co- 
efficients Pn(L) along the iso-energetic trajectory. An 
exception is the case when the energy scale of all the pro- 
cesses involved is such that only transitions between the 
ground state (n = 1) and the first excited state (n = 2) 
are possible. In this effective two-level spectrum, com- 
bining Eq. ([15p with the normalization condition Eq.©, 
the trajectory for the iso-energetic process is described 
by the following relation 



Pi(L) 



E 2 (L 1 )-E 2 (L) g!(Li)-£; 2 (Li) 
E 1 (L)-E 2 (L) E 1 {L)-E 2 {L) 



(19) 



with p 2 (L) = 1 —pi(L) after the normalization condition 
Eq. © • The heat exchanged by the system with the envi- 
ronment during an iso-enegetic trajectory connecting the 
initial and final states L\ L 2 , for the case of a Dirac 
particle, is given by the expression 

-Q x _> a - [E?(I*) + (^f(Lx) - E?(L X )) 
■ E?(L 2 )-E?(L 2 y 
_E?(L 1 )-E?(L 1 )_ 
■Eg(L 2 ) mc 2 4 
■Eg (Li)mc 2 -f 



xpi(Li)] In 
+(mc 2 ) In 



mc 



mc 2 - 



■ 2(mc 2 ) In 



-E?{L 2 ) 
-E?{L l )_ 



'hi 



(20) 



where E%(L) was defined in Eq.©. 

For the effective two-level system previously described, 
we shall conceive a cycle, as depicted in Figure 1, which 
starts in the ground state with p\{La) = 1. Then, 
the system experiences an iso-energetic expansion from 
La — > Lb > La- Then, it experiences an iso-entropic 
expansion from Lb — > Lc > Lb, then an iso-energetic 
compression L c — > L D < L c , and finally it goes back to 
the initial ground state through an iso-entropic compres- 
sion Lb, — > La- 

We shall assume that the final state after the iso- 
energetic process La — > Lb corresponds to maximal ex- 
pansion, that is, the system ends completely localized 
in the excited state n = 2. In this condition, Eq. (fT9"|) 
reduces to 



px{L B )=Q, P2 {L B ) = 1. 



(21) 



The condition of total energy conservation between the 
initial and final states connected through an iso-energetic 
process, for maximal expansion, leads to the equation 



p 1 (L A )E 1 (L A ) = p 2 {Lb)E 2 [Lb) 



(22) 



where pi(L A ) — P2{Lb) = 1 for maximal expansion. 
Therefore Eg. (|2"2")l . given the Dirac spectrum Eq.©, im- 
plies that Lb /La = 2. 



The heat released to the environment along this first 
stage of the cycle is calculated after Eq. (j2"(3} , 



Qa^b = E°(L A )\n 



E 2 D (L A )-E?(L J 



Eg{2L A )-Eg{2L / 



+ mc 2 In 



'2 
mc 



2 1 Ef(L, 



E i D (2L A ) 



(23) 



The next process is an iso-entropic expansion, charac- 
terized by the condition p 2 (Lb) = pi(Lc) = E We shall 
define the expansion parameter a = Lc/Lb > 1. The 
work performed during this stage, with Lb = 2L A (as 
discussed before), is calculated from Eq. (fT4|) . 




(24) 



The cycle continues with a maximal compression pro- 
cess from Lc = 2uLa to Lb = uLa under iso-energetic 
conditions. The energy conservation condition is in this 
case similar to Eq. (|22|) . with p 2 {Lc) = Pi(Ld) = E 
The heat exchanged by the system with the environment 
along this process, applying Eq. (|2"0)) , is given by the ex- 
pression 



Q 



C^D 



Eg (2aL A )ln 



E 2 L> ( a L A )-E 1 D (aL A ) 



+mc In 



Eg(2aL A )-E 1 D (2aL A ) 
nc 2 + Eg{aL A ) \ 



E°(2aL A ) J 



(25) 



where E®(L) was defined in Eq.©. The last path along 
the cycle is an adiabatic process, which returns the sys- 
tem to its initial ground state withpi(L£>) = px(L A ) = 1. 
The work performed during this final stage, as obtained 
by applying Eq. (|14[) . is given by 



W D - 




(26) 



It is interesting to check that the work along the two iso- 
entropic trajectories cancels, that is Wb^,c + Wd^,a = 0. 
Therefore, the efficiency of the cycle is defined by the 
ratio 



V D = 1 _ VC^D 



Q 



A-¥B 



(27) 



When substituting the corresponding expressions from 
Eq.([2"3")) and Eq.([2"5|) into Ea.([2T]). we obtain the explicit 
analytical expression 



,D 



= 1 - 



In 


1 1+9(qL a /2) f 


6(aL A /2) 


-6{aL A ) \ v 


0" 


4 


l+8(2aL A ) \ 


9(aL A )~ 


0{2aL A ) ] 




In 


1 l+0(L A /2) 


(6(L A /2) 


-0(L A )\ e( - LA y 




4 1+8(2L A ) 


1 8(L A )- 


6(2L A ) j 





(28) 
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FIG. 4: (Color online)Pictorial description of the Force versus 
width of the potential well, for the Quantum Carnot Cycle 
described in the text. 



Here. 



have defined the function 
with A = 27rft/(mc) the 
It is important 



6{L) = 
Compton 
that in the 



y/1 + (A/2i)2 

wavelength. It is important to remark tnat m 
non-relativistic limit X/L — > 0, the expression in Eq.([27l 
reduces to the known Schodinger limit 



lim n 

X/L-yO 



1 - 1/a 2 



(29) 



In the " ultra- relativistic" case of a massless Dirac 
fermion, X/L — > oo, the efficiency converges to the length 
independent limit 



lim 

A/L^c 



1 - 1/a. 



(30) 



The trend of the efficiency is shown in between both lim- 
its in Figure 3. It is worth to remark that, since the 
expansion parameter a — Lq/Lb > 1, the efficiency in 
the strict non-relativistic (Schrodingcr) regime X/L — > 
Eq. ([29"]) is the highest possible one. This is also clear from 
the asymptotics of the curves displayed in Fig. 3, where 
the " ultra-relativistic" limit corresponding to massless 
fermions (X/L — > oo) indeed represents the less efficient 
regime for a fixed expansion parameter a. 



V. THE QUANTUM CARNOT CYCLE 

In this section, we shall discuss the quantum version 
of the Carnot cycle, as applied to the statistical ensem- 
ble of Dirac singlc-fermion systems under consideration. 
The thermodynamic cycle which defines the correspond- 
ing heat engine is composed of four stages or trajectories 
in Hilbert's space: Two isothermal and two iso-entropic 
processes. 

In the first stage, the system is brought into contact 
with a thermal reservoir at temperature Tu- By keeping 
isothermal conditions, the width of the potential well is 
expanded from La — > Lb- Since thermal equilibrium 
with the reservoir is assumed along this process, the von 



Neumann entropy of the system achieves a maximum for 
the Boltzmann distribution [l6|, [TtJ 

Pn (L, fa) = [Z(L, faT 1 e-^W, (31) 

with /3 = (k B T)- 1 , and the normalization factor is given 
by the partition function 



Z(L,/3) = J2e 



X 



'K^pmc 2 ). (32) 



Here, the second expression, as shown in Appendix, is 
the continuum approximation to the discrete sum, valid 
in the physically relevant regime A<L. Here, Ki(x) is 
a Besscl function of the second kind. 

From a similar analysis as in the previous section, we 
conclude that the heat exchanged by the system to the 
thermal reservoir is given by 



Q 



.4-s-B 



>ln'^ 



dp n (L,fa) 



dL 



dL 



Z(L A ,p H ) 



dfa 



H JZ(L B ,fa) 
\Z{L A ,fa) 



" 2 lnf ^ 



(33) 



In the second line, we have done integration by parts, 
and we made direct use of the definition Eq. (|30"]) of the 
partition function. The final result follows from substi- 
tuting the explicit expression for the partition function 
Eq.©. 

Similarly, during the third stage of the cycle, the sys- 
tem is again brought into contact with a thermal reser- 
voir, but at a lower temperature Tq < Tu- Therefore, 
the probability distribution of states in the ensemble is 
p n {L,f3c), as defined in Eq . (|31"jl . but with Tc instead of 
Th- The heat released to the reservoir during this stage 
is given by the expression 



Qc^d = mc In 



(34) 



The second and fourth stages of the cycle consti- 
tute iso-entropic trajectories. In order to analyze these 
stages, we shall derive the "equation of state" for the 
statistical ensemble of single-particle systems. When 
substituting the Boltzmann distribution p n (/3,L) = 
\Z(f3,L)}~ ex-p(—f3E®(L)) into the expression for the 
von Neumann entropy Eq. (|10p . we obtain the relation 

(35) 



S = T~ 1 E + k B \nZ(P,L). 



Here, E = (H) is the ensemble-average energy, as de- 
fined by Eq. (jTTJ) . The equation of state is obtained from 
Eq.® as 



F = - (§) s = kBT ^ nmL) 
k R T 



L 



(36) 
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In the last line, we have used the explicit analytical ex- 
pression Eq.(f3"2"j) for the partition function to calculate 
the derivative. The equation of state Eq. ([36| reflects that 
the ensemble of systems behaves as a one-dimensional 
ideal gas. This is not surprising, since the ensemble- 
average energy is given by 



E = (H) = -—]nZ{p,L) 

= mc 2 1-l-^lnK^z 
\ dz 



(37) 



where we have substituted explicitly the expression for 
the partition function, and in the final step we defined 
z = mc 2 . Here, Ki{x) is a Bcsscl function of the second 
kind. Eq. (|37)) shows that the ensemble average energy of 
the system is a function of the temperature solely, from 
which the ideal gas equation of state follows as a natural 
consequence. We can thus define the "specific heat" at 
constant length, which after Eq. ([37|) is given by 



Cr 



dE 
df 



§ = k BZ ^K l{z) , (38) 



where z = j3mc 2 . It is interesting to remark that, 
based on the asymptotic behaviour of the Bessel func- 
tions K n (z) ~ v /7r/2z~ 1 / 2 exp(— z) + 0{z~ 1 ), the spe- 
cific heat defined in Eq. (|38)l presents the asymptotic limit 
Cl — > k B /2 when fc^T <C mc 2 . This is the well known 
result for a classic non-relativistic ideal gas in one di- 
mension. This feature and the general temperature de- 
pendence of the ensemble specific heat is displayed in 
Fig. 6. The change in the ensemble averaged energy of 
the system, for a general process, is dE = TdS — FdL. 
Since the ensemble-average energy is a function of tem- 
perature only, the differential equation for an iso-entropic 
trajectory (dS = 0) is(l9j 



dE = C L dT = -FdL = -k R T^. 



Separating variables, after some algebra we obtain 



d 2 

z-j-^ In Ki(z) 



dL 
T' 



(39) 



(40) 



Integrating Eq. (|4Tfl) between initial conditions (z Q ,L ) 
and final conditions (z, L), we have 



L 



Ki(z ) 
Ki(z) 



(41) 



Here, we made use of the Bessel function identity 
K 1 (z) = Kq(z) — z~ x Kx(z), with z = (3mc 2 . It is inter- 
esting to check that in the non-relativistic limit z » 1, 
given the asym ptotic behaviour of the Bessel functions 
K n (z) ~ N /^72z- 1 / 2 exp(-z) + 0{z- x ), Eq.gU) reduces 
to iT -1 / 2 = cat. for the iso-entropic trajectory. On the 
other hand, in the " ultra- rclativistic" limit of a massless 




FIG. 5: (Color online)Pictorial description of the Force versus 
width of the potential well for the Quantum Carnot Cycle. 



Dirac particle with z — > 0, the iso-entropic trajectory is 
given by the equation LT^ 1 = ant. 

We are now in conditions to discuss the second and 
fourth stages of the Carnot cycle. The second stage of 
the process corresponds to an iso-entropic trajectory, pa- 
rameterized in differential and integral form by Eqs. (|40[) 
and (|4"Tj) . respectively. The work performed by the sys- 
tem during this process is given by 



W B - 









FdL = 


Jl b 




= —mc 2 


r c d 2 



La 



dL 
k B T— 



L 



In K\(z)dz 



(42) 



Here, in the second line we have used the differential 
equation defining the iso-entropic trajectory, Eq.(|4Ul). 
Evaluating the integral in Eq. (|4"2"|) . we explicitly obtain 



WB^c = k B (Tc-T H ) 



Kp(zc) _ K (z H ) 



(43) 



The fourth and final stage of the cycle also corresponds 
to an iso-entropic trajectory Ld — > La, and the work per- 
formed by the system against the external applied force 
is obtained similarly as in Eq. (|4"2")l , 



W D - 



FdL 



k B {Th - T c 



(44) 



K (z H ) K Q (z c ) 



K^zh) K x {z c ) 



Clearly, after Eqs.([43]l and (|35|). we have W B ^c + 
Wd^a — 0, and hence the contribution of the work along 
the iso-entropic trajectories vanishes. 

From the equation for the iso-entropic trajectory 
Eq. (l4Tj) , we conclude that the length ratios are deter- 
mined by the temperatures of the thermal reservoirs, 



L C _L D _ K 1 (z H ) cZC $$g. 
L B L A K 1 (z c ) 



(45) 



From Eq. (|4"5"|) . we also obtain La/L b = L B /Lc- Sub- 
stituting this relation in the expression for the efficiency 
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FIG. 6: (Color online) Specific heat of the statistical ensemble 
of single Dirac-fermion systems, after Eq.(j38}, as a function 
of temperature. 



of the cycle, from Eq. and Eq. we have 
_ l _ Qc^d _ l _ T c ^.(Ld/Lc) 



Q 



T H HLb/La) 



= i-*£. 



(46) 



Therefore, we have recovered the expression for the effi- 
ciency identical to the classical Carnot cycle. 



VI. CONCLUSIONS 

By considering as a "working substance" the statisti- 
cal ensemble for a Dirac single-fermion system trapped in 
a one-dimensional potential well, we have analyzed two 
different schemes for a quantum heat engine. The first, 
that we have referred to as the E-Carnot cycle, consists 
of two iso-entropic and two iso-energetic trajectories. We 
obtained an explicit expression for the efficiency of this 
cycle and showed that our analytical result, in the non- 
relativistic limit X/L — > 0, reduces to the corresponding 
one for a Schrodinger particle, as reported in the litera- 
ture Our results also indicate that the efficiency for 
this E-Carnot cycle is higher in the non-relativistic region 
of parameters, that is for L>A, when comparing at the 
same compressibility ratios a = Lc/Lb > 1. An excep- 
tion is the case of massless Dirac fermions, with A = oo, 
where it is not possible to achieve non-relativistic condi- 
tions. This is of potential practical interest for graphene 
systems, where conduction electrons are indeed described 
as massless chiral Dirac fermions [TlTfTij . 

As a second candidate for a quantum heat engine, we 
discussed a version of the Carnot cycle, composed by 
two iso-thermal and two iso-entropic trajectories. In or- 
der to achieve iso-thermal conditions, we consider that 
the single-fermion system is in thermal equilibrium with 
macroscopic reservoirs at temperatures Tq < Tjj, respec- 



tively. Therefore, the statistical ensemble under consid- 
eration is described by the density matrix p = er^ H /Z, 
with Z\p,L] = Tre~ p6 the partition function. We 
showed that the statistical properties of the ensemble are 
such that an equation of state can be defined, as well as a 
specific heat, in analogy with a classical ideal gas in one 
dimension. We obtained the equation for the iso-entropic 
trajectory, which in the non-relativistic limit ksT <C mc 2 
reduces to the classical result LT~ X / 2 = ant. On the 
other hand, we also showed that in the " ultra-relativistic" 
limit of a massless Dirac fermion, as for instance conduc- 
tion electrons in graphene, the iso-entropic trajectory is 
denned by the equation LT^ 1 = cat. We also showed 
that the efficiency for the Quantum Carnot cycle satis- 
fies the same relation that the classical one in terms of the 
temperatures of the thermostats, that is rj = 1 — Tc/Tjj. 
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Appendix 

We here discuss the continuum approximation to the 
discrete partition function Eq.(j30[), 



71=0 



Here, let us define the discrete variable x n = n\/(2L). 
The spacing between two consecutive values is given by 
Ax = x n +i — x n = X/(2L), and hence the expression for 
the partition function can be written as 

Z (L, p) = ^f^Ax e~ 0mc2 (V^- 1 ) 

71 = 

For physically relevant sizes of the potential well, we shall 
have X/L <C 1, which allow us to take the continuum limit 
in the sense of a Riemann sum for the previous equation, 

* Jo 

= 2 j± r dte -^c 2 (co S h ( t)-i) cosh(t) 

A Jo 



2L 
T 



Here, in the second line we have made the substitution 
x = sinh(t), and Ki(x) is a Bessel function of the second 
kind. 
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